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(Dated: November 10, 2015) 

Using 2.92 fb -1 of electron-positron annihilation data collected at a center-of-mass energy of 
y/s = 3.773 GeV with the BESIII detector, we present an improved measurement of the branching 
fraction B(D + —> uje + v e ) = (1.63 ±0.11 ±0.08) x 10 -3 . The parameters defining the corresponding 
hadronic form factor ratios at zero momentum transfer are determined for the first time; we measure 
them to be ry = 1.24±0.09±0.06 and r 2 = 1.06±0.15±0.05. The first and second uncertainties are 
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statistical and systematic, respectively. We also search for the decay D + —» (j>e + v e . An improved 
upper limit B(D + —» (j>e + Ve.) < 1.3 x IGF 5 is set at 90% confidence level. 

PACS numbers: 13.20.Fc,14.40.Lb 


Charm semileptonic decays have been studied in detail 
because they provide essential inputs of the magnitudes 
of the Cabibbo-Kobayashi-Maskawa (CKM) elements 
\V c d\ and \V„\ [1, 2], and a stringent test of the strong 
interaction effects in the decay amplitude. These effects 
of the strong interaction in the hadronic current are 
parameterized by form factors that are calculable, for 
example, by lattice QCD and QCD sum rules. The 
couplings \V CS \ and \V cc l\ are tightly constrained by the 
unitarity of the CKM matrix. Therefore, measurements 
of charm semileptonic decay rates and form factors rig¬ 
orously test theoretical predictions. Both high-statistics 
and rare modes should be studied for a comprehensive 
understanding of charm semileptonic decays. 

For D —y Vtv transitions (where V refers to a vector 
meson), the form factors have been studied in the decays 

D + —> K*°e + v e [3] and D + — > p°e + u e [4]. The decay 
D + —> u>e + v e was first observed by the CLEO-c experi¬ 
ment, while the corresponding form factors have not yet 
been measured due to limited statistics [4], The decay 
D + —> Loe + v e can proceed through the tree-level diagram 
shown in Fig. 1. Its transition rate depends on the 
charm-to-down-quark coupling \V c d\, which is precisely 
known from unitarity of the CKM matrix. Neglecting 
the lepton mass, three dominant form factors contribute 
to the decay rate: two axial (Ai,A 2 ) and one vector ( V ) 
form factor, which are functions of the square of the 
invariant mass of the lepton-neutrino system q 2 . 

The decay D + — > (pe + v e has not yet been observed. 
The most recent experimental search was performed by 
the CLEO Collaboration in 2011 with a sample of an in¬ 
tegrated luminosity of 818 pb _1 collected at the 03770) 
resonance. The upper limit of the decay rate was set to 
be 9.0 x 10 -5 at the 90% confidence level (C.L.) [5]. Since 
the valence quarks ss of the <f> meson are distinct from 
those of the D meson ( cd ), this process cannot occur in 
the absence of u>-<f> mixing or a non-perturbative “weak 
annihilation” (WA) contribution as shown in Fig. 2 [6, 7]. 
A measurement of the branching fraction can discrimi¬ 
nate which process is dominant. For example, a study of 
the ratio of Df —> we + zz e and Df — > 4>e + v e [6] concludes 
that any value of B(Df uje + v e ) exceeding 2 x 10 -4 
is unlikely to be attributed to u}-(f> mixing, and would 
provide evidence for non-perturbative WA effects [7]. A 
search for the decay D + — > <f>e + v e is helpful, since its 
dynamics is similar to that of the decay Df —> we + ^ e . 

In this paper, we report an improved measurement of 
B{D + —> uje + v e ) and the first form factor measurement 
in this decay. Furthermore, an improved upper limit for 
B(D + (j)e + v e ) is determined. Charge conjugate states 
are implied throughout this paper. Those decays are 
studied using a data sample collected with the BESIII 
detector which corresponds to an integrated luminosity 



d _«_ d 

FIG. 1. Feynman diagram representing the tree-level charged 
current process D + —> uoe + u e . 



FIG. 2. Feynman diagram representing the WA process 
D + —¥ (j>e + v e . 


of 2.92 fb 1 at the •0(3770) resonance [9]. 

The BESIII detector is a spectrometer located at 
BEPCII, which is a double-ring e + e _ collider working 
at the center-of-mass energy range from 2 to 4.6 GeV. 
The cylindrical core of the BESIII detector consists of a 
helium-based multi-layer drift chamber (MDC), a plas¬ 
tic scintillator time-of-flight system (TOF), and a Csl 
(Tl) electromagnetic calorimeter (EMC), which are all 
enclosed in a superconducting solenoid magnet providing 
a 1.0 T magnetic field. The solenoid is supported by an 
octagonal flux-return yoke with modules of resistive plate 
muon counters interleaved with steel. The momentum 
resolution for charged particles at 1 GeV/c is 0.5%, 
and the resolution of the ionization energy loss per unit 
path-length ( dE/dx ) is 6%. The EMC measures photon 
energies with a resolution of 2.5% (5%) at 1 GeV in the 
barrel (end cap). The time resolution of the TOF is 80 
ps in the barrel and 110 ps in the end cap. A detailed 
description of the BESIII detector is provided in Ref. [10]. 

The tagging technique for the branching fraction mea¬ 
surements of semileptonic decays was first employed by 
the Mark-Ill collaboration [11] and later applied in the 
studies by CLEO-c [4, 12]. The presence of a D + D~ 
pair in an event allows a tag sample to be defined in 
which a D~ is reconstructed in one of the following six 
hadronic decay modes: -K + 7 t - 7t - , K, K°n~, 
Kgir 7r°, KgTT + TT n , and K + K 7r .A sub-sample is 
then defined in which a positron and a set of hadrons 
are required recoiling against the tag D meson, as a 
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signature of a semileptonic decay. The yields of tag 
and signal are expressed as 7V t l ag = 2 N D jjBl e[ and 
N si g = 2N d£>B ta g #si<4g, s i> respectively, where N dB is 
the total number of DD pairs produced, i indicates a 
tag mode, V t l ag is the number of observed tag events in 
mode i, AA g is the number of semileptonic candidates, 
B\ ag is the branching fraction of tag mode i. B s \ is the 
branching fraction of the semileptonic decay, e[ ag is the 
reconstruction efficiency of a tag mode, and e[ ag sl is the 
reconstruction efficiency of the semileptonic decay with 
a tag mode. Thus, B s \ can be expressed as 


H^tagetag.slAtag ’ 
i 

where N s ; g is the total signal yield in all six tag modes. 

Charged tracks are reconstructed using hit informa¬ 
tion from the MDC. The tracks are required to satisfy 
| cos < 0.93, where 6 is the polar angle with respect 
to the beam axis. Tracks (except for Kg daughters) are 
required to originate from the interaction point (IP), i.e. 
their point of closest approach to the interaction point 
is required to be ±10 cm along the beam direction and 
1 cm transverse to the beam direction. Charged particle 
identification (PID) is accomplished by combining the 
dE/dx and TOF information to form a likelihood 
(* = e/n/K) for each particle hypothesis. A K± (tt^) 
candidate is required to satisfy Ck > C K ( C v > Ck)- 
For electrons, we require the track candidate to satisfy 
c +c‘+c K > as well as E/p £ [0.8,1.2], where E/p 
is the ratio of the energy deposited in the EMC to the 
momentum of the track measured in the MDC. In order 
to take into account the effect of final state radiation and 
bremsstrahlung, the energy of neutral clusters within 5° 
of the initial electron direction is assigned to the electron 
track. The Kg candidates are reconstructed from pairs 
of oppositely charged tracks, which are assumed to be 
pions and required to have an invariant mass in the range 
m n + n - £ [0.487,0.511] GeV/c 2 . For each pair of tracks, 
a vertex-constrained fit is performed to ensure that they 
come from a common vertex. 

To identify photon candidates, showers must have min¬ 
imum energies of 25 MeV in the barrel region (| cos0| < 
0.80) or 50 MeV in the end cap region (0.86 < j cos0| < 
0.92). To exclude showers from charged particles, a 
photon candidate must be separated by at least 20 ° from 
any charged track with respect to the IP. A requirement 
on the EMC timing suppresses electronic noise and en¬ 
ergy deposits unrelated to the event. The 7 r° candidates 
are reconstructed from pairs of photon candidates by 
requiring the invariant di-photon mass to fulfill m 77 £ 
[0.115,0.150] GeV/c 2 . Candidates with both photons 
coming from the end cap region are rejected due to poor 
resolution. 

The D~ tag candidates are selected based on two 
variables: A E = Ed — Abeam, the difference between 
the energy of the D~ tag candidate (£x>) and the 


TABLE I. Tag yields of the D six hadronic modes and their 
statistical uncertainties. 


Tag mode 

Nt a K 

D~ K + tv~'it~ 

809425 ± 906 

D~ — >■ A' + 7T _ 7T _ 7r 0 

242406 ± 599 

D~ —>■ K% n~ 

100149 ± 321 

D~ -> K° s n~Tr° 

226734 ± 575 

D~ —>■ Kg n + Tv~n~ 

132683 ± 489 

D~ —>■ K+K~n~ 

70530 ± 325 

Total 

1581927 ±1399 


beam energy (Abeam), and the beam-constrained mass 
M bc = y/ A 2 eam /c 4 - |pd| 2 /c 2 , where p D is the measured 
momentum of the D~ candidate. In each event, we 
accept at most one candidate per tag mode per charge, 
and the candidate with the smallest |AA| is chosen. The 
yield of each tag mode is obtained from fits to the M^ c 
distributions [13]. The data sample comprises about 
1.6 x 10 6 reconstructed charged tag candidates (Table I). 

Once a D~ tag candidate is identified, we search for 
an e + candidate and an w —> n + Tr~TT° candidate or a 
(j> —> K + K~ candidate recoiling against the tag. If 
there are multiple ui candidates in an event, only one 
combination is chosen based on the proximity of the 
7 r + 7 r _ 7 T° invariant mass to the nominal ui mass [14]. 
The invariant mass m n + 7r - n o £ [0.700,0.840] GeV/c 2 
and m K + K ~ £ [1.005,1.040] GeV/c 2 are required for 
w and (j) candidates, which correspond to three times 
of the u {(/) mass resolution (±3cr), respectively. For 
the decay D + —> we + i/ e , backgrounds arise mostly from 
D + — y K*°e + v e , K*° —y Ag7r 0 , I\g —> n + n~ process, 
and the invariant mass of the charged pions is required 
to be outside the aforementioned Kg mass region. This 

: | : Q 

requirement rejects about 70% of the D + —> K e + v e 
background events. 

After tag and semileptonic candidates have been com¬ 
bined, all charged tracks in an event must be accounted 
for. The total energy of additional photon candidates, 
besides those used in the tag and semileptonic candi¬ 
dates, is required to be less than 0.250 GeV. Semileptonic 
decays are identified using the variable U = A m j ss — 
c|Pmiss|, where A miss and p miss are the missing energy 
and momentum corresponding to the undetected neu¬ 
trino from the D + meson semileptonic decay, which are 
calculated by A m j ss = Abeam A^,( 7 ^ A e , p m iss = 
— {pt&g + PuiU) + Ve) in the center-of-mass frame, where 
A u (0) (A e ) and p^u) (p e ) ar e the energy and momentum 
of the hadron (electron) candidate. To obtain a better U 
resolution, the momentum of the tag D~ candidate ptag 
is calculated by p t ag = Ftag[(A be am/c) 2 - M|,c 2 ] 1/2 [ 15 ], 
where ptag is the unit vector in the direction of the tag 
D~ momentum, and Mp is the world average value of D 
meson mass [14]. The correctly reconstructed semilep¬ 
tonic candidates are expected to peak around zero in the 
U distribution. A GEANT4-based [16] Monte Carlo (MC) 
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TABLE II. Tag efficiencies (etag) and signal efficiencies in¬ 
cluding a tag (etag.ai) in percent, as determined by the MC 
simulation, and their statistical uncertainties. 


Tag 

; mode 


etag 


£tag,sl (^) 

etag 

,sl 

(</>) 

D~ 

K+ n~ 

TV 

51.07 

± 

0.02 

11.22 + 0.10 

9.04 

± 

0.09 

D~ 

->• K + TT~ 

'n~n° 

25.13 

± 

0.02 

5.15 + 0.09 

4.38 

± 

0.08 

D~ 

~^K° S 7T- 


54.40 

± 

0.05 

11.70 + 0.32 

9.69 

± 

0.29 

D~ 

—>• K% 7T 

7T° 

29.24 

± 

0.02 

6.13 + 0.11 

5.34 

± 

0.10 

D~ 

->• K° s n+ 

7V~7V~ 

37.61 

± 

0.04 

7.28 + 0.18 

5.96 

± 

0.16 

D~ 

-+ K+K 

~TV~ 

41.12 

± 

0.06 

8.97 + 0.29 

7.63 

± 

0.27 


simulation is employed, and events are generated with 
KKMC+EVTGEN [17, 18] to determine the efficiencies in 
Eq. (1), as shown in Table II. All selection criteria and 
signal regions are defined using simulated events only. 

The yield of the decay D + —> oje + v e is obtained from 
a fit to the U distribution combining all tag modes, as 
shown in Fig. 3. The signal shape is described by the 
shape from the signal MC simulation convoluted with 
a Gaussian function whose width is left free in the fit 
to describe the resolution difference between MC and 
data. The background model consists of two components: 
peaking and non-peaking backgrounds. Peaking back- 

ground arises mostly from the decay D + —> K e + u e , 
K —> Kg-rr 0 , Kg —f 7 r + 7 r _ ; its U distribution is modeled 
with MC simulation. The largest contribution to the non¬ 
peaking backgrounds is from the DD process, while the 
remaining background events are from the non -DD, qq , 
t + t - , initial state radiation 7 J/i^ and r yip(2S) processes. 
The non-peaking component is modeled with a smooth 
shape obtained from MC simulations. In the fit to data, 
the yield of the peaking background is fixed to the MC 
expectation, while that of the non-peaking background 
is left free in the fit. The signal yield is determined by 
the fit to be JV s i g = 491 ± 32. The absolute branching 
fraction of the decay D + -+ ue + v e as listed in Table III 
is obtained using Eq. (1). 

The U distribution for the decay D + —> (j)e + v e with all 
tag modes combined is shown in Fig. 4. The signal region 
is defined as [—0.05,0.07] GeV, which covers more than 
97% of all signal events according to MC simulations. 
No significant excess of signal events is observed, and 
there are only 2 events in the signal region. A simulation 
study indicates that the backgrounds arise mostly from 
D + — > (j> 7T + 7r° and D + — > <fiir + processes. The number of 
background events is estimated to be 4.2 ± 1.5 via large 
statistics MC samples. The upper limit is calculated 
by using a frequentist method with unbounded profile 
likelihood treatment of systematic uncertainties, which 
is implemented by a C++ class TROLKE in the ROOT 
framework[19]. The number of the observed events is 
assumed to follow a Poisson distribution, and the number 
of background events and the efficiency are assumed to 
follow Gaussian distributions. The resulting upper limit 
on B(D + —► (j)e + v e ) at 90% C.L. is determined to be 



FIG. 3. Fit (solid line) to the U distribution in data (points 
with error bars) for the semileptonic decay D + — » ioe + u e . The 
total background contribution is shown by the filled curve, 
while the peaking component is shown by the cross-hatched 
curve. 



U(GeV) 


FIG. 4. The U distribution for the semileptonic decay 
D + —» (j>e + v e in data (points with error bars) and signal MC 
simulation with arbitrary normalization (solid histograms). 
The arrows show the signal region. 


1.3 x 10 5 , as shown in Table III. 

With the double tag technique, the branching fraction 
measurements are insensitive to systematics from the tag 
side since these are mostly cancelled. For the signal side, 
the following sources of systematic uncertainty are taken 
into account, as summarized in Table IV. The uncer¬ 
tainties of tracking and PID efficiencies are well 

studied by double tagging DD hadronic decay events. 
The uncertainties in eF tracking and PID efficiency are 
estimated with radiative Bhabha events. The uncertainty 
due to the 7 r° reconstruction efficiency is estimated with 
a control sample D° —y K~n + TT 0 by the missing mass 
technique. The uncertainty due to imperfect knowledge 
of the semileptonic form factors is estimated by varying 
the form factors in the MC simulation according to the 
uncertainties on the measured form factor ratios in the 
decay D + -+ oje + u e as discussed below. For the decay 
D + -+ (j)e + v e , the signal MC produces phase-space dis¬ 
tributed events, and therefore uses a constant form fac¬ 
tor. To evaluate the corresponding systematics, the form 
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TABLE III. Measured branching fractions in this paper and 
a comparison to the previous measurements [4, 5]. For the 
D+ —» uje + v e decay, the first uncertainty is statistical and 
the second systematic. 

Mode This work Previous 

we+u e (1.63 ± 0.11 ± 0.08) x 10~ 3 (1.82 ± 0.18 ± 0.07) x 10~ 3 
cpe+ve < 1-3 x 10~ 5 (90%C.L.) < 9.0 x 10~ 5 (90%C.L.) 


factor is varied by a reweighting technique [8] . The world 
average values of B(uj —>■ 7r + 7r~7r° ) and B(cj) —► K + K~) 
are (89.2±0.7)% and (48.9±0.5)%, respectively, and their 
uncertainties are assigned as systematic uncertainties due 
to the input branching fractions in the MC simulation. 
The limited MC statistics also leads to a systematic 
uncertainty. The uncertainties associated with the ui 
and 4> mass requirements are estimated using the control 
samples D° —> ujK~tt + and D + —> cj>ir + , respectively. 
The Kg rejection leads to an uncertainty on the signal 
efficiency of the decay D + we + v e , which is studied by 
the control sample D° —> ujK~tt + . The uncertainty due 
to the extra shower veto is studied with double hadronic 
tags. For the decay D + —> (f>e + v e , the uncertainty due 
to the signal window requirement as shown in Fig. 4 is 
estimated by the control sample D + —> K e + u e , K —> 


K~tt + . In the fit to the U distribution in the D + —> 
we + v e decay, the uncertainty due to the parametrisation 
of the signal shape is estimated by varying the signal 
shape to a Crystal Ball function [20]. The uncertainty 
due to the fit range is estimated by varying the fit range. 
The uncertainty due to the non-peaking background is 
estimated by modeling this component with a third-order 
Chebychev function, and the uncertainty associated with 
the fixed peaking background normalization is estimated 
by varying it within its expected uncertainty. All of those 
estimates are added in quadrature to obtain the total 
systematic uncertainties on the branching fractions. 

The differential decay rate of D + —> uje + v e can be ex¬ 
pressed in the following variables as illustrated in Fig. 5: 
to 2 , the mass square of the 7 r 7 T 7 r system; g 2 , the mass 
square of the ev e system; 9 i, the w helicity angle [ 21 ], 
which is the angle between the uj decay plane normal (n) 
in the 7 T 7 r 7 r rest frame and the direction of flight of the u> 
in the D rest frame; 02, the helicity angle of e, which is 
the angle between the charged lepton three-momentum 
in the ev e rest frame and the direction of flight of the 
ez'e system in the D rest frame; x, the angle between the 
decay planes of those two systems. 

For the differential partial decay width, only the P- 
wave component is taken into consideration and the 
formalism expressed in terms of three helicity amplitudes 
H+{q 2 ), H_(q 2 ), and H 0 (q 2 ) is [4, 22, 23]: 


dr 

dq 2 d cos 9 \d cos 92d\dm- K ^- K 


= 8(Lr G2F]Vcdl2l j^ B{u} 7r7r7r )i^ w ( m ^)i 2 

[(1 + cos 6 » 2 ) 2 sin 2 9 1 \H + (q 2 ,m^)\ 2 

+ (1 - cos 9 2 ) 2 sin 2 9i\H-(q 2 , m 7r7r7I -)| 2 +4sin 2 9 2 cos 2 9i\H 0 (q 2 , m^T,)] 2 
+ 4sind 2 (1 + cos 0 2 ) sin#! cos#i cos \H+{q 2 , m n7Tn )H Q (q 2 , 

— 4sin# 2(1 — cos # 2 ) sin#i cos#i cos \H-{q 2 , TO 7 r 7 r 7 r )iJo(9 2 I w W7r7r ) 

- 2 sin 2 # 2 sin 2 #1 cos 2 \H+ (q 2 , (q 2 , to^)]. 


( 2 ) 


where Gf is the Fermi constant, is the momentum of 
the to in the D rest frame, B(u> —> 7 T 7 r 7 r) is the branching 
fraction of ui —> mnr, is the invariant mass of the 

three pions, and 6 W(to wtw ) is the Breit-Wigner function 
that describes the w line shape. The helicity amplitudes 
can in turn be related to the two axial-vector form factors 
4 li. 2 (# 2 ) and the vector form factor V(q 2 ): 


H± (q 2 )=MA 1 (q 2 )T2^V(q 2 ) 


M 

t2 


H 0 (q 2 )=- - — m„„ — q 2 )MAi{q 2 ) (3) 

2 to^VS 

,M > 2 


-4- 


M 


- Mq 2 )] 


where M = Mo + . For the q 2 dependence, a single 


pole parameterization [24] is applied: 


V(q 2 ) 


^( 0 ) 

1 - q 2 /m^ 


AiAq 2 ) 


^1,2(0) 
l-q 2 /m 2 A ’ 


( 4 ) 


where the pole masses my and toa are expected to 
be close to = 2.01 GeV/c 2 and = 

2.42 GeV/c 2 [14] for the vector and axial form factors, 
respectively. The ratios of these form factors, evaluated 
at q 2 = 0 , ry = and r 2 = , are measured in 

this paper. 

According to the fit procedure introduced in Ref. [3], 
a five-dimensional maximum likelihood fit is performed 
in the space of m 2 , g 2 , cos#i, cos #2 and %• The signal 
probability density function is modeled with the phase- 
space signal MC events reweighted with the decay rate 
(Eq. 2) in an iterative procedure. The experimental 
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TABLE IV. Summary of systematic uncertainties on the 
branching fraction measurements. 


Source 

B(D+ -a weU e ) B{D+ ->■ <t>e+v e ) 

Tracking 

3.0% 

3.0% 

K/ 7T PID 

1 . 0 % 

1 . 0 % 

e PID 

3.2% 

3.4% 

7r° reconstruction 

1 . 0 % 

- 

Model of form factor 

1 . 0 % 

1 . 2 % 

decay rate 

0 . 8 % 

1 . 0 % 

MC statistics 

0.7% 

0.9% 

cmass window 

0.9% 

0.4% 

Kg veto 

0 . 2 % 

- 

Extra shower veto 

0 . 1 % 

0 . 1 % 

Signal window 

- 

0.4% 

Fit range 

0.4% 

- 

Signal shape 

0 . 6 % 

- 

Peaking background 

0 . 8 % 

- 

Non-peaking background 

0.4% 

- 

Total 

5.1% 

5.0% 



FIG. 5. Definitions of the helicity angles in the decay D + —> 
ljW + , u) —» 7 r + 7 r _ 7 r°, W + —> e + v e for the three-body ( 8 1 ) 
and two-body ( 82 ) D + -daughter decays, where both angles 
are defined in the rest-frame of the decaying meson. 


acceptance is taken in consideration using this technique. 
Large signal MC samples are generated to reduce the 
systematic uncertainty associated with the MC statistics. 
The background is modeled with the MC simulation 
and its normalization is fixed to the expectation. Using 
simulated events with known ry and rqi, we verify that 
this procedure can reliably determine the form factor 
ratios. Figure 6 shows the to 2 , q 2 , cos 6 h, cos 62 and \ 
projections from the final fit to data. The fit determines 
the form factor ratios to be ry = 1-24 ± 0.09 and 
r 2 = 1.06 ±0.15. 

For the form factor measurement in the decay D + —>- 
coe + u e , the following sources of systematic uncertainties 
are taken into account, as summarized in Table V: The 
uncertainty associated with the unknown q 2 dependence 
of the form factors is estimated by introducing a double 
pole parameterization [25]. The uncertainty due to the 


TABLE V. Summary of the absolute systematic uncertainties 
in the form factor measurement of the decay D + —> oje + v e . 


Sources 

rv 

r 2 

q 2 dependence 

0.05 

0.03 

Background model 

0.02 

0.02 

Pole mass assumption 

0.01 

negligible 

Fitting shift 

0.02 

0.02 

Total 

0.06 

0.05 


background model is estimated by varying the back¬ 
ground normalization with its statistical uncertainty. No 
events from the non-resonant decay D + —> tt + n~ n°e + v e 
are observed, the influence of this decay on the form 
factor therefore can be neglected. To estimate the un¬ 
certainty associated with the pole mass assumption, we 
vary the pole mass my by ±100 MeV/c 2 . A small 
shift is observed with the presence of background, and 
this is treated as possible bias in the form factor fit¬ 
ting procedure. Adding all systematic uncertainties in 
quadrature, the form factor ratios are determined to be 
r v = 1.24 ± 0.09 ± 0.06 and r 2 = 1.06 ± 0.15 ± 0.05, 
respectively. 

In summary, using 2.92 fb _1 of e + e _ annihilation data 
collected at i/s = 3.773 GeV with the BESIII detector, 
we have measured the form factor ratios in the decay 
D + —> uje + v e at q 2 = 0 for the first time: ry = Ai(o) = 
1.24 ± 0.09 ± 0.06, r 2 = = 1.06 ± 0.15 ± 0.05, 

and determined the branching fraction to be B(D + —> 
ioe + v e ) = (1.63 ± 0.11 ± 0.08) x 10 -3 , where the first and 
the second uncertainties are statistical and systematic, 
respectively. This is the most precise measurement to 
date. We have also searched for the rare decay D + —> 
(f>e + v e and observe no significant signal. We set an upper 
limit of B(D + —>• (j)e + v e ) < 1.3 x 10 ~ 5 at the 90% C.L., 
which improves the upper limit previously obtained by 
the CLEO Collaboration [5] by a factor of about 7. 
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FIG. 6. Projections of the data set (points with error bars), the fit results (solid histograms) and the sum of the background 
distributions (filled histogram curves) onto (a) m 2 , (b) q 2 , (c) cos#i, (d) cos 62 and (e) X". 
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